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Two-nucleon systems are shown to exhibit large scattering lengths in strong magnetic fields at 
unphysical quark masses, and the trends toward the physical values indicate that such features 
may exist in nature. Lattice QCD calculations of the energies of one and two nucleons systems 
are performed at pion masses of ~ 450 and 806 MeV in uniform, time-independent magnetic 
fields of strength |B| ~ 10 1!) 10 2tl Gauss to determine the response of these hadronic systems 

to large magnetic fields. Fields of this strength may exist inside magnetars and in peripheral 
relativistic heavy ion collisions, and the unitary behavior at large scattering lengths may have 
important consequences for these systems. 

PACS numbers: 11.15.Ha, 12.38.Gc, 13.40.Gp 


In most physical situations, external electromagnetic 
(EM) fields have only small effects on hadronic and nu¬ 
clear systems, whose structure and dynamics are dom¬ 
inated by the internal strong interactions arising from 
Quantum Chromodynamics (QCD) and internal EM in¬ 
teractions. However, there are specific situations involv¬ 
ing extremely large EM fields, created either naturally 
in astrophysical environments or in particle colliders, for 
which the effects of external fields are important. In mag¬ 
netars, high magnetic field rotating neutron stars [T] , sur¬ 
face magnetic fields are observed up to 0(1O 14 ) Gauss 
(for reviews, see e.g. Ref. (21 0 ), and it is conjectured 
that interior magnetic fields reach up to 0(1O 19 ) Gauss 
[3]. In heavy ion collisions, the currents produced by rel¬ 
ativistic nuclei lead to large magnetic fields within the 
projectiles, particularly during (ultra-)peripheral colli¬ 
sions [5]- F i s estimated that fields of C>(10 19 ) Gauss 
are experienced by the nuclei during the femtoseconds 
of the nuclear crossings 5j. Neither of these environ¬ 
ments are easy to probe in a controlled way, and the 
detailed behavior of nuclei in such fields is an open ques¬ 
tion. As a step toward exploring nuclei in these extreme 
magnetic environments, we present the results of calcu¬ 
lations of the effects of uniform, time-independent mag¬ 
netic fields on two-nucleon (as well as two-hyperon) sys¬ 


tems performed with the underlying quark and gluon de¬ 
grees of freedom. We find that such fields can poten¬ 
tially unbind the deuteron and significantly modify the 
nucleon-nucleon (NN) interactions in the [ Sq channel. At 
the unphysical quark masses where the calculations are 
performed, the scattering lengths in both the 3 S\- 3 Di 
and ',S'o channels diverge at particular values of the field 
strength. Near these values, the low energy dynamics 
of these systems will become unitary. The trends seen 
towards the physical values of the quark masses suggest 
that this feature may exist in nature in some of these sys¬ 
tems. The prospect of such resonant behavior in nuclear 
systems is exciting and it will be important to incorpo¬ 
rate this effect into models of magnetars and heavy ion 
collisions in which the relevant field strengths are probed. 

Before presenting the results of our calculations, it is 
interesting to consider phenomenological expectations for 
the behavior of such systems. 1 For small, constant mag¬ 
netic fields, the responses of the two-nucleon systems be- 


1 Significant effort has been devoted to understanding the nature 
of the QCD vacuum in strong magnetic fields (see Ref. [9] for 
a review), but effects specific to hadronic systems are not well 
studied. 
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yond their charges are governed by their magnetic mo¬ 
ments if the system has spin, and otherwise by their mag¬ 
netic polarizabilities. The deuteron has a magnetic mo¬ 
ment such that in a magnetic field in the 2 direction 
the j z = +1 component is positively shifted in energy 
with respect to the breakup threshold 2 and so an ap¬ 
proach toward unbinding in a magnetic field is plausible. 
However, higher order responses to the magnetic field 
may be important, and at intermediate field strengths, 
|eB| ~ m 2 , significant deviations from linearity should be 
anticipated. In the opposite limit of extremely large mag¬ 
netic fields, where |eB| Aq CD , the asymptotic freedom 
of QCD implies [10] that the eigenstates evolve towards 
weakly-interacting up and down quarks in Landau lev¬ 
els. Hence, as the magnetic field tends to infinity, the 
ground states of dilute systems tend to threshold. When 
the density of the system is also large and comparable to 
the scale of Landau orbits, more exotic phases may occur 
(see Ref. |Hj for a review). 

In this work, the numerical technique of Lattice QCD 
(LQCD) is applied to study two-nucleon systems in uni¬ 
form, time-independent background magnetic fields, fol¬ 
lowing methods used in previous studies of the magnetic 
moments [12] and polarizabilities m of nucleons and 
light nuclei up to atomic number A = 4. To understand 
the phenomenological effects of the strong fields in nu¬ 
clear environments, a first task is to ascertain the ef¬ 
fects on the two-nucleon interactions. Two particle scat¬ 
tering phase shifts can be accessed in LQCD from the 
volume dependence of two-nucleon energies (the Luscher 
method PH US]), but here a simpler approach is under¬ 
taken in which only the bound states of the two-nucleon 
sector are addressed. 3 The primary goal of these calcu¬ 
lations is to investigate how the binding energies of the 
two-nucleon states respond to applied magnetic fields. 

LQCD calculations were performed using two ensem¬ 
bles of gauge-field configurations generated with a clover- 
improved fermion action [Tfi] and the Luscher-Weisz 
gauge action Cl- The first ensemble had Nf = 3 de¬ 
generate light-quark flavors with masses tuned to the 
physical strange quark mass, producing a pion of mass 
m n ~ 806 MeV, and used a volume of L 3 xT = 32 s x 48. 
The second ensemble used Nf = 2 + 1 quark flavors 
with the same strange quark mass and degenerate up and 
down quarks with masses corresponding to a pion mass of 


2 While the deuteron magnetic moment is positive, it is less than 
the sum of the neutron and proton magnetic moments. In a 
potential model the difference is due to the d-state admixture 
into the predominantly s-wave deuteron wave function, while in 
NN effective field theories (EFTs) this is encapsulated in short- 
distance two-nucleon interactions with the magnetic field. 

3 At unphysically large values of the light quark masses, both the 
deuteron and dineutron are bound, as are various two baryon 
hypernuclei ns- 


m„ ~ 450 MeV and a volume of L 3 x T = 32 3 x 96. Both 
ensembles had a gauge coupling of /3 = 6.1, corresponding 
to a lattice spacing of a ~ 0.11 fm. The ensembles con¬ 
sisted of ~ 1, 000 gauge-field configurations at the SU(3) 
point and ~ 650 configurations at the lighter pion mass, 
each taken at intervals of 10 hybrid Monte-Carlo trajec¬ 
tories. We have extensively studied these ensembles in 
previous works, and have found that the finite-volume 
effects to both the single nucleon and two-nucleon bound 
state energies are small mm- 

As in Refs. [T21 Hoi [2U] . background EM {Uq( 1)) gauge 
fields were implemented through the gauge-links, 
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that give rise to uniform magnetic fields along the in¬ 
direction. These were multiplied onto each QCD gauge 
field in each ensemble (separately for each quark flavor of 
charge Q q ). The combined QCD+EM gauge fields were 
used to calculate up-, down-, and strange-quark propa¬ 
gators, which were then contracted to form the requi¬ 
site nuclear correlation functions using the techniques 
of Ref. [2T]. To ensure periodicity, n € Z, and the 
values n = 0,1,—2,3,4,—6,12 were used on the SU(3) 
symmetric ensemble, while n = 0,1, —2,4 were used on 
the m,r ~ 450 MeV ensemble. The corresponding field 
strengths are quantized as |eB| = 67r|fi|/(aL) 2 , giving 
a field of C>(10 19 ) Gauss for n = 1. On each con¬ 
figuration, quark propagators were generated from 48 
uniformly distributed Gaussian-smeared sources for each 
magnetic field. For further details of the production at 
the SU(3)-symmetric point, see Refs. 12] [18] H2] and in 
particular, Ref. [13]. Analogous methods were used for 
the light mass ensemble. 

This work focuses on the dineutron, the diproton, and 
the maximal \j z \ = j = 1 spin state of the deuteron, 
all of which remain isolated, sub-threshold states in the 
presence of a magnetic field. The I z = j z = 0 neutron- 
proton systems with (j = 1;J = 0) and (j = 0;/ = 1) 
mix in a magnetic field and have been considered previ¬ 
ously in Ref. [2U] to determine the cross section for the 
radiative capture process np —> dr/. States with the quan¬ 
tum numbers of h = n,p,nn, pp 1 du i—i are accessed 
from correlation functions Ch(t; B) = (0|x?i(t)x?i(0)|0)B 
computed in the presence of the background magnetic 
field B from source and sink interpolating operators with 
the requisite quantum numbers, as discussed in detail in 
Ref- [13]. Representative correlation functions for the 
heavier mass ensemble can be found in Ref. [13] for each 
hadron/nucleus and background magnetic field. Ratios 
of these correlation functions to those without the mag¬ 
netic field, Rh{t; B) = Ch(t\ B)/C/,(i; 0), are also shown 
in Ref. (13] . and are used to extract the magnetic mo¬ 
ments and polarizabilities of the respective systems. For 
the m n ~ 450 MeV ensemble, the ratios behave in a 
qualitatively similar manner and the signals are of com¬ 
parable quality. As the central focus of this study is on 
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the difference between the effect of the field on the two- 
nucleon systems and on the nucleons in isolation, the 
further ratios 

6R A (t;B) = / H R h (t;B) , (2) 

' heA 

are of primary importance. In this expression, A refers 
to the composite system and the product is over its con¬ 
stituent nucleon correlator ratios (e.g., for A = dj z =+1 
the contributions are from p 4 and n4). The late time ex¬ 
ponential decay of this ratio is dictated by the binding 
energy of the system in the presence of the field na. 
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Fig- 0 shows these ratios for the m^ ~ 450 MeV en¬ 
semble along with the results of single exponential fits to 
time ranges in which the individual correlation functions 
entering the ratios are consistent with single exponen¬ 
tial behavior. As discussed in Ref mi, multiple different 
interpolating operators are investigated for each state in 
this study and the resulting differences are used to gauge, 
in part, the systematic uncertainty. In the figures below, 
we focus on a particular set of interpolating operators for 
clarity but have verified that other choices of interpola¬ 
tors provide consistent results. The analogous results for 
the heavier mass ensemble are presented in Ref. m ■ 
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FIG. 1: Correlator ratios defined in Eq. for the nn, 
the j z = +1 deuteron and pp systems for field strengths 
n = 1, —2, 4, for the m n ~ 450 MeV ensemble. The bands cor¬ 
respond to the exponential fit and its statistical uncertainties 
associated with the shown fit interval. Systematic uncertain¬ 
ties from the choice of fit range are separately assessed. 


The energy shifts 

A^(n) = SE A ( B) - Y, * E h( B) (4) 

heA 

in the dineutron and deuteron ( j z = +1) channels are 
shown in Figs. [2] and [3j respectively. As the strength of 
the applied magnetic field is increased, the ground state 
energies of the systems are shifted closer to threshold, 
and at a given field strength it appears that the states 
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FIG. 2: Response of the binding of the dineutron system to 
applied magnetic fields. The upper panel shows the result at 
ro, = 806 MeV, while the lower panel is for m n = 450 MeV. 
The shaded regions correspond to the envelopes of successful 
fits to the energy shifts using linear and quadratic polynomials 
in n 2 to data points in the corresponding range indicated by 
the shaded region. The horizontal bands indicate the binding 
threshold. 



FIG. 3: Response of the binding of the j z = +1 state of 
the deuteron to applied magnetic fields. The shaded regions 
correspond to the envelopes of successful fits to the energy 
shifts using polynomials in h of up to A th (2 nd ) order for 
the m-K = 806 (450) MeV ensemble. The horizontal bands 
indicate the binding threshold. 


unbind. For the deuteron, this behavior is not clearly 
resolved at the lighter mass because of the uncertainties. 
The approach to threshold and subsequent turnover is 
seen at both quark masses in the dineutron system, and 
the point of minimum binding decreases as the quark 
mass is lowered, ~ 6 at ~ 806 MeV and 

fxn™ ax ' 1 ~ 3 at ~ 450 MeV. The dineutron is unbound 
in nature and the present results suggest that magnetic 
effects would push the system further into the contin¬ 
uum. On the other hand, it is possible that the deuteron 
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FIG. 4: Response of the binding of the diproton to ap¬ 
plied magnetic fields. The shaded regions correspond to the 
envelopes of successful fits to the energy shifts using polyno¬ 
mials in h of up to 4 th (2 nd ) order for the m n = 806 (450) 
MeV ensemble. The horizontal bands indicate the binding 
threshold. 



FIG. 5: The energy splitting in the H dibaryon channel at 
ra 7 r = 806 MeV. The horizontal band indicates the binding 
threshold. 

could be unbound by the presence of magnetic fields of 
strength comparable to those expected in magnetars and 
heavy ion collisions, potentially modifying the dynamics 
of those systems. A particularly interesting aspect of the 
behavior in both of these channels is the approach to the 
unitary regime in which the binding energies decrease 
to zero and consequently the scattering lengths diverge. 4 
In atomic physics, such behavior is routinely used to in¬ 
vestigate the universal physics that emerges in systems 
interacting near unitarity [22] , but they have not been 
observed in nuclear physics. 

The energy shifts of the diproton are shown in Fig. [4] 
For this system, the extracted energies are not as cleanly 
determined as for the dineutron, but a trend toward 
strengthening attraction is seen at both quark masses 


4 It is expected that the range of the interaction (set by hadronic 
scales) is only weakly affected by the magnetic field, so the vol¬ 
ume effects in the two-nucleon systems are not expected to be 
unmanageable even as the scattering length diverges. 


as the field strength increases. This is interesting in light 
of a recent suggestion J23j that the diproton can over¬ 
come the Coulomb repulsion and form a bound state in 
sufficiently large magnetic fields. A naive extrapolation 
of the slope of the shift linearly in to); indicates that for 
a field of |eB| ~ 10 17 Gauss, corresponding to n ~ 0.01, 
the additional attraction is enough to bind the diproton 
system. While such a result would be interesting, fur¬ 
ther calculations at lighter quark masses are necessary 
to refine the extrapolation. 

Two baryon systems containing strange quarks have 
also been investigated. Figure 5 shows the energy split¬ 
tings of the ground state in the channel with the quan¬ 
tum numbers of two A-baryons, which contains a deeply 
bound FZ-dibaryon at heavier quark masses [25125]. This 
channel exhibits a slight reduction of the binding energy 
for intermediate field strengths, comparable in size to 
that of the dineutron system, but does not exhibit res¬ 
onant behavior in the range of field strengths that are 
probed as the binding energy is significantly larger. 

Discussion: Having found significant changes in the 
binding of two-nucleon systems immersed in strong mag¬ 
netic fields at two values of unphysical quark masses, it 
is conceivable that similar modifications occur in nature. 
To solidify this discussion, calculations would need to 
be performed at or near the physical quark masses and 
the continuum and infinite volume limits would require 
careful investigation. 5 While the responses of these sys¬ 
tems can as yet only be estimated at the physical quark 
masses, the calculated trends provide an interesting start¬ 
ing point to consider possible consequences. To this end, 
it is conjectured that two-nucleon systems will exhibit 
unitary behavior, with the deuteron unbinding in a large 
magnetic field and the diproton system becoming bound. 
On the other hand, the dineutron will be pushed further 
into the continuum as the field strength increases. In¬ 
terestingly, it may be possible to find values of the field 
strength and quark masses where all NN states are at 
threshold simultaneously, realizing the low energy confor¬ 
mal symmetry postulated by Braaten and Hammer [26j . 
Given the observed behavior of bound states, it is nat¬ 
ural to expect that the NN scattering phases shifts and 
mixing angles will also be modified at a similar level in 
such fields. These modifications would be interesting to 
probe in future LQCD calculations utilizing the Liischer 
method mm to analyze the spectra of NN systems. 

In ultra-peripheral heavy ion collisions, one can spec¬ 
ulate that the reduced binding between pairs of nucle¬ 
ons, along with the reduction in the nucleon mass, will 
increase the size of each nucleus as they interact with 


5 Based on studies of binding energies on these and other related 
ensembles, we are confident that the current calculations do not 
suffer from large volume or scaling artifacts. 
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the field of the other nucleus. Ignoring other potential 
effects, purely geometrical considerations will result in 
larger than expected interaction cross-sections that will 
increase with the collision energy for a given impact pa¬ 
rameter and potentially larger fluctuations in collision 
cross sections. However, considering the transient nature 
of such a collision, and the difference between the internal 
time-scales associated with a rearrangement of the nu¬ 
cleons comprising each nucleus and that of the collision, 
more detailed analyses must be performed before even a 
qualitative understanding can be established. The effects 
of large magnetic fields in magnetars through the mag¬ 
netic moments of nucleons and electrons have been con¬ 
sidered through a number of model approaches [a emi], 
and in some cases lead to significant modifications. The 
more complicated effects from magnetic shifts in binding 
and hadronic interactions likely also induce significant 
modifications that deserve further investigation. 
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